Selective coherent destruction of tunneling in a quantum-dot array 
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The coherent manipulation of quantum states is one of the main tasks required in quantum 
computation. In this paper we demonstrate that it is possible to control coherently the electronic 
position of a particle in a quantum-dot array. By tuning an external ac electric field we can selectively 
suppress the tunneling between dots, trapping the particle in a determined region of the array. 
The problem is treated non-perturbatively by a time-dependent Hamiltonian in the effective mass 
approximation and using Floquet theory. We find that the quasienergy spectrum exhibits crossings 
at certain field intensities that result in the selective suppression of tunneling. 
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The search for a solid-state based quantum computer 
device has attracted a lot of interest in the physics com- 
munity recently. The ability to manipulate a quantum 
state and measure it is one of the most important and 
pressing challenges to be addressed in real implementa- 
tions. Fortunately great advances in this area have been 
achieved recently. A good example of that are the Rabi 
oscillations observed in exciton states of self-assembled 
quantum dots (QDs). 1 ! 2 ! 3 ! 4 

Our work here is based on the coherent destruction 
of tunneling (CDT) of a driven two-level system^*S*i in 
which the tunneling of one particle in a symmetric dou- 
ble well potential is suppressed for some special value 
of frequency and field intensity. We propose a system in 
which one can selectively suppress the tunneling between 
individual quantum dots tuned by the external ac electric 
field. This effect provides one with a different experimen- 
tal handle to control a quantum mechanical system. By 
suitable variation of the frequency and applied ac field 
amplitude, one can precise the location of one electron in 
the multidot array. To demonstrate this effect, our model 
employs an effective mass nearest-neighbor tight-biding 
approximation (NNTB)&2ii£ The dynamics is analyzed 
using Floquet theory and the direct integration of the 
time-dependent Schrddinger equation. 

The Hamiltonian for an electron in an array of identical 
QDs under a strong ac field within NNTB is written as 

H = T e (aj +1 aj + h.c.) + eFdjdjdj cos(cj£ + <fi), 

3 3 

(1) 

where T e is the hopping matrix element, Oj and aj are, 
respectively, the electron creation and annihilation op- 
erator in the dot j, e is the electronic charge, F is the 
field intensity, d is the separation between dots, w is the 
field frequency, and is the phase of the drive field. This 
phase, first thought to be a non-relevant factor, is in fact 
quite an important parameter in the CDT at smaller fre- 
quency, as we have recently showniAi Indeed, a phase 
cj> = 7r/2 produces a much better dynamic localization 



than any other possible phase cj>, and consequently it is 
our choice here. 

Since H is periodic in time (H(i) = H(t + r), 
where r = 2tt/uj is the period) we can make use of 
the standard Floquet theory^ 2 , and write the solutions 
of the time-dependent Schrodinger equation as ip(t) = 
exp(—iet/K)u(t), where u(t), the so-called Floquet state, 
is also periodic in time with the same period r, and e 
is the Floquet characteristic exponent or quasienergy, 
which can be obtained from the eigenvalue equation 



(2) 



H - ih— ) u(t) = eu(t). 



Note that this equation is similar to the time- 
independent Schrodinger equation with TL = H — ifidt 
playing the role of a time-independent Hamiltonian. Us- 
ing this analogy we can explore a combined dynamic par- 
ity operation: z — > —z;t —> t + r/2, under which the 
operator TL is invariant. As a consequence, each Flo- 
quet state is either even or odd under this operation; 12 
Quasienergies of different "dynamic parity" may cross, 
while an avoided crossing is expected as a function of ex- 
ternal parameters, such as the field intensity, for states 
with the same parity. 

In a symmetric double QD the levels present strictly 
different parity, and they may exhibit crossings with field 
unless the symmetry is broken. If we vary the field inten- 
sity, and the quasienergy levels cross, a CDT is expected, 
since a splitting between levels comes from this inter-dot 
tunneling. This CDT result is well known in the high 
frequency limit, and it occurs at field values satisfying 
the zeroes of the Bessel function, J (eFd/huj) — 0. As 
we have shown in Ref. 0, and extended by Crefheld,- 13 
CDT may also occur at lower frequency with decreasing 
degree of localization and at different intensities of the 
driven field. The dynamics of the system can be dras- 
tically different, depending on the numbers of dots. We 
note that for this effect to be observed, the frequency 
of the driven field cannot be excessively high, compared 
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FIG. 1: (Color online) First Brillouin zone of quasienergies 
(in units of Ttw) as function of ac field intensity (eFd/hui), 
for T e /huj — 0.2. Inset shows amplification of collapse region, 
showing crossings and anticrossing that follow a definite pat- 
tern according to dynamical symmetries. The two vertical 
lines indicate field intensity Fi and F2 of pair crossings and 
are used later in numerical simulation. 



with the tunneling probability between nearest dots, as 
otherwise the level splitting will be so small as to be un- 
noticeable. It should also be pointed out that in the limit 
of an infinite number of identical QDs the levels form a 
miniband with bandwidth of 4T e , and since there is an 
infinite number of levels in that range, they arc infinitcsi- 
mally separated, so that even at low frequencies the levels 
collapse at the zeroes of the Bessel function. Notice that 
the applicability of the NNTB in that limit system has 
some limitationsii^*i^*i& The system considered here has 
just a few QDs and we assume weak coupling between 
them, so that the NNTB can be applied with confidence. 
Notice also that we are considering a different kind of 
dynamic localization since we do not focus only in the 
probability to find the particle in the same initial state, 
but in a determined region of the dot array that may 
not be necessarily the initial localized state. In fact the 
particle stays confined to a determined multidot region 
of the array. 

In our analysis, we focus on the dynamics of four QDs. 
However, we show that similar results can be obtained 
for other numbers of QDs as well. In Fig. ^ we show 
the quasienergy spectrum as a function of the field inten- 
sity for the ratio T e /Tiw = 0.2 exhibiting the strong de- 
pendence on field amplitude and the expected collapse of 
the spectrum near field values for which Jo(eFd/huj) = 0. 
Notice however, as shown in the inset, that the region of 
the level collapse is in fact a region with several crossings 
and anticrossings following a well-defined pattern given 
by the dynamic parity. The two vertical lines in the in- 
set are shown to indicate the two field values, F1 and 
F2 , where pairs of levels cross and are the choices of field 
intensity used for the numerical simulations below. 

Our numerical simulation is done by direct integration 
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FIG. 2: (Color online) Time evolution of probability to find 
one particle in each one of the QDs in a four-dot array for 
eF\d/Tiuj ~ 2.38, and starting the system with the particle 
in: (a) dot 1, and (b) dot 2. Lower panel is a schematic 
representation of the dynamics of the system, as seen in the 
respective time evolution of the occupation probability of each 
dot. The full circle represents the position of the particle 
at time zero, crosses indicate the suppression of tunneling 
through that barrier and arrows indicate that tunneling is 
possible. 



of the time-dependent Schrodinger equation and followed 
by calculations of the occupation probability for differ- 
ent choices of initial conditions. In Fig. [2Ja) we show 
the time evolution of the system for the field intensity 
F\ corresponding to the first vertical line in the inset of 
Fig- (eFid/hu! ~ 2.38), assuming that at time zero the 
particle is in the first QD. Notice that the particle is basi- 
cally frozen in that dot, which means that for this choice 
of field we can effectively suppress the tunneling in the 
barrier between dots 1 and 2. This is what the bottom 
cartoon in Fig.[3represents, where the cross in the barrier 
indicates that the tunneling is not allowed and that the 
effective interdot tunneling amplitude T e —> 0. On the 
other hand, if we start the system with the particle in dot 
2, Fig. Gib) shows that the particle can tunnel back and 
forth between dots 2 and 3, while the effective tunneling 
is not only suppressed between dots I and 2, but between 
3 and 4 also. The bottom panel represents the dynamics 
for this initial condition. We can then conclude that for 
this first choice of field intensity one can selectively block 
the tunneling of the electron through the outer barriers, 
so that the only open barrier for tunneling is the one in 
the middle 

Let us now explore the second choice of field inten- 
sity F2, the second vertical line in the inset of Fig. ^ 
(eF2d/huj ~ 2.4). This small but discernible tuning of 
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FIG. 3: (Color online) Time evolution of the probability to 
find a particle in one of the QDs for eFd/hw ~ 2.40, starting 
the system with the particle in: (a) dot 1, and (b) dot 2. 
Lower panel is a schematic representation of the dynamics. 
Here, the middle barrier tunneling is suppressed. 



the field intensity yields a completely different result. At 
that point there are two crossings in the spectrum. In 
Fig. Ota) , we show the results when we start the system 
with the particle in dot 1, and in Fig. IHIb) starting in 
dot 2. The picture for both initial conditions does not 
change at all, except for a 7r-phase in the oscillation (that 
is provided by the initial state) . The effective tunnclings 
can be summarized in the respective bottom panels of 
that figure, and we can see that now we can block the 
tunneling between dots 2 and 3. Similar behavior and 
conclusions are reached for initial conditions in either dot 
3 or 4. 

These examples represent an exciting result since by 
simply tuning the field intensity from F± to F2, one can 
choose which barrier is allowed or blocked for tunneling. 
Notice that this dot or site selectivity is achieved despite 
the ac field being applied to the entire structure, which 
suggests interesting applications. Notice, moreover, that 
the amount of field intensity needed to change this condi- 
tion is relative to the frequency of the driving field, i.e., 
the ratio T e /huj. It is also interesting that the separa- 
tion between the two tunneling suppressing fields is not 
a monotonic function of T e /huj. Figure ^a) shows how 
these two crossing points change with the ratio T e /hu). 
The solid line shows the first crossing at field intensity F\ , 
and the dashed line is the result for F^. Notice that in the 
high frequency regime (smaller ratio T e /hu)) the crossings 
occur very close to each other and at field intensities sat- 
isfying the condition Jo(eFd/huj) = 0. As a result, all 
barriers present a suppression of tunneling at the same 
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FIG. 4: (Color online) (a) Crossing in the collapse region as 
indicated by vertical lines in the inset of Fig. Q as function 
of the ratio T e /huj. In the high frequency limit (lower ratio 
T e /huj) the crossings occur basically at the same point (zeroes 
of Bessel function) . Inset shows amplification of quasienergy 
spectrum for T e /huj = 0.7. (b) Escape probability for a par- 
ticle initially in one region of the quantum-dot array. Solid 
red line is for the equivalent first crossing and represents the 
escape probability of the QD 1, considering this dot as the ini- 
tial condition. Dotted blue line continuation shows F\ when 
the crossing disappears. Dashed green line is for the second 
crossing and represents escape probability of the first two QDs 
assuming that the system starts in QD 1. 



value of field intensity in the high frequency regime 
For lower frequencies and given T e , the crossings become 
separated and different barriers are selectively closed at 
different fields. At even lower frequency T e /Tiuj > 0.4, the 
first crossing disappears (indicated by the dotted blue 
line continuation) as the coupling between quasienergy 
levels one and tree, and two and four increases. This re- 
sults in an effective gap in the quasienergy spectrum (see 
inset of Fig.QJa) for T e /Tioj — 0.7) that basically kills the 
CDT for that value of field intensity and frequency. 

To better understand how good is this CDT at values 
of field given by F± and F2 we can monitor the maximum 
probability for the particle tunneling out of the region of 
dots we are considering (see Ref. 12). Lower values for 
this represent a good CDT. Figure^b) shows the escape 
probability for field intensities F\ and F2 provided by the 
crossings in Fig. 12 a), as function of the ratio T e /hu) (fre- 
quency). The solid line is the result for the first crossing 
Fi , and is the escape probability out of dot 1 , defined as 
the maximum value reached by P2 + P3 + Pa, assuming 
that the system starts with the particle in dot 1. The 
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FIG. 5: (Color online) Collapse region of the first Brillouin 
zone of quasienergies (in units of hui) as function of the ac 
field intensity (eFd/hui), for T e /hu — 0.2 for: (a) an array 
of six QDs and (b) an array of five QDs. Lower cartoons are 
the schematic representation of the dynamic of the respective 
arrays of QDs at field intensity at crossing points given by F\ 
and F'2. 

same result is clearly obtained if the system starts with 
the particle in dot 2 and one monitors the maximum 
probability to find it in dot 1. The dotted blue line that 
follows the solid red line in Fig.QJb) is the result for field 
values F\ at frequencies when there are no more cross- 
ing, but a gap in the quasienergy spectrum. Notice this 
goes to unity quickly as the effective tunneling becomes 
possible as a result of having a gap in the quasienergy 
spectrum. The dashed green line in Fig. HJb) i s the re- 
sult for the second level crossing field F2, representing 
the escape probability of the dots 1 and 2 (since for this 
choice of field we suppress the tunneling between dots 2 
and 3), which is the maximum value reached by P3 + P4 
starting the system with the particle in either dots 1 or 



2. Notice that for frequencies around T e /hw ~ 0.3 the 
system still presents a good CDT for both crossings, and 
they occurs at field intensity with significant separation. 

Our analysis was given for four QDs, but the same be- 
havior can be obtained for other finite number of dots. 
The quasienergy spectrum presents two crossings, F\ and 
F2, in the collapse region as we can see in Fig.[SJa) for six 
and (b) for five QDs, respectively. The first crossing F\ 
basically suppresses the tunneling between dots 1 and 2, 
and between the two last QDs. This result is schemati- 
cally represented in the bottom cartoons for F\ for either 
five or six QDs. The second crossing in the spectrum, F2, 
suppress the tunneling between pairs of quantum dots, 
for example, in the six QDs cases, the tunneling between 
dots 2 and 3, and 4 and 5 as schematically represented in 
the lower bottom panel of Fig.|5[a). An interesting result 
appears for odd numbers of QDs, since in this case there 
is a different number of even and odd states. This fact 
allows for the dot in the middle to become a true trap 
for the particle, since tunneling in both directions can be 
suppressed. For example, the case of five QDs for the 
second crossing, F%, as schematically represented by the 
lower bottom cartoon of Fig.^d), suppress the tunneling 
between the dots 2 and 3, and 3 and 4. So, if we could 
start the system, or manipulate it, in a state localized in 
that dot and choose this value of field intensity (F2) we 
could in principle trap the particle in the middle of the 
quantum-dot array. 

In conclusion, we have shown that it is possible to sup- 
press selectively the tunneling between quantum dots by 
a simple tuning from F\ to F2, the intensity of an applied 
ac field. With this tool one could manipulate the posi- 
tion of the particle in a quantum-dot array, and assist in 
its initialization and control. 
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